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A quasi  one-dimensional  foil  bearing  jnodel  for  high  speed  rotating  magnetic  heads  has  been 
formulated.  Side  flow  is  •‘beared"  as  a leakage  factor.  Steady  state  foil  inertia  effects  are 
included.  The  effects  of  head  protrusion  side  leakage  factor  as  well  as  speed  parameter  are 
discussed.  Of  particular  interest  is  the  effect  of  rotational  speed  in  relation  to  the  propagation 
speed  of  transverse  waves  in  the  foil. 
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characteristic  dimension  in  side  flow  direction 
bending  stiffness  of  foil 

dimensionless  deviation  of  foil  contour  — 1 

local  fluid  film  thickness 

asymptotic  film  thickness  away  from  the  head 

dimensionless  film  thickness  h/h* 

local  pressure  under  foil 

ambient  pressure 

pressure  under  foil  in  vicinity  of  head 

drum  radius 

time 

foil  tension  per  unit  width 
foil  velocity  (in  +x  direction) 
drum  velocity 

dimensionless  foil  contour  w/h* 

foil  contour  (function  of  x and  t) 

head  protrusion 

bump  contour  (function  of  x) 

circumferential  position  along  drum 

side  flow  parameter,  proportional  to  (6/a)**2 

head  width  at  base 

speed  parameter  (Eq.  12) 

air  viscosity 

foil  mass  per  unit  area 

dimensionless  time  (Eq.  10) 

dimensionless  circumferential  position  x/6 
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I INTRODUCTION 

1 

I High  density  tape  recorders  in  computer  and  instrumentation  appli- 

< cations  utilize  a geometry  schematically  shown  in  Figure  1.  A drum  in  which  a 

magnetic  head  is  imbedded,  rotates  at  a high  speed.  The  tape  is  wrapped  around 
the  drum  and  is  floated  over  it  with  a fairly  thick  clearance  (order  of  25  microns). 
Due  to  the  rotational  speed,  the  tape  tends  to  float  also  at  a considerable  separa- 
tion from  the  head.  It  is  desired,  however,  that  the  head  interface  the  tape  at  a 
close  proximity  (a  fraction  of  a micron).  This  task  is  achievable  by  especially 
contouring  and  grooving  the  head  to  promote  reduction  of  the  hydrodynamic 
film  thickness  to  the  desired  level.  The  effects  of  side  flow  are  thus  of  signifi- 
cance. Due  to  the  high  speed  of  drum  rotation,  the  effects  of  foil  inertia  are 
also  of  importance.  In  fact,  one  may  think  of  the  head  as  forming  a concentrated 
travelling  load. 

The  objective  of  the  present  report  is  to  approximate  on  a quasi 
one-dimensional  basis  as  many  effects  as  possible  of  the  configuration  described 
in  Figures  1-3.  Specifically,  the  effects  of  head  contour  and  protrusion,  side  flow 
(represented  by  a "smeared"  leakage  factor),  foil  inertia,  tension  and  speed  will 
be  discussed.  Results  as  to  the  effect  of  the  various  parameters  will  be  reported. 

A selection  from  the  relevant  literature  is  mentioned  below.  The 
elastic  effects  of  travelling  concentrated  disturbances  in  plates  or  foils  have  been 
studied  by  Kenney  [1],  Reismann  [2],  Bogy,  Greenberg  and  Talke  [3],  and 
Adams  and  Bogy  [4].  Langlois  [5]  studied  a quasi  two  dimensional  membrane 
model  with  lubrication  effects  included.  Dais  and  Barnum  [6]  looked  into  the 
f effect  of  a rotating  drum  with  a flat.  Greenberg  [7]  investigated  the  interface 

of  a magnetic  head  with  a flexible  disk.  Eshel  examined  the  effect  of  a travelling 
external  pressure  feed  hole  [8].  Stahl,  White  and  Deckert  [9]  published  work  on 
dynamic  effects,  taking  foil  inertia  into  consideration.  In  the  experimental  area, 
the  main  available  work  is  by  Albrecht,  Laenen  and  Chua  Lin  [10]  who  reported 

» - 

measurements  of  foil  contours  in  response  to  travelling  bumps. 
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ANALYSIS 

Figure  1 describes  schematically  the  model  being  studied.  A drum 
of  radius  r^,  with  a bump  of  a specified  contour  w^(x)  rotates  with  a surface 
velocity  Here,  x is  a circumferential  coordinate  fixed  in  the  drum.  The 
foil  wraps  the  drum,  and  its  contour  Wf(x,t)  as  well  as  the  foil  to  drum  separ- 
ation h(x,t),  are  the  subject  of  the  investigation.  An  essentially  planar  model 
is  being  studied  here.  It  is  assumed,  however,  that  part  of  the  circumferential 
flow  is  removed  by  side  flow  into  appropriate  grooves.  An  approximate  math- 
ematical model  for  the  side  flow  is  hypothesized  later.  The  force  balance  across 
the  foil  is: 

'■»  ' ax'* 

^ ’’'^ax'a4-  <r  ' 

The  one-dimensional  Reynolds  equation  in  these  coordinates  with  a local  leakage 
correction  factor  included  is: 

where 

SIDE  LEAKAGE  3 ^ ^ 1(1  P 
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Figure  1 Schematic  View  of  a Drum,  Bump  and  Foil 

TYPICAL  SIDE  FLOW  CONTROLLING  DIMENSION 


Figure  2 View  of  Lateral  Contour  of  a Grooved  Bump 
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In  the  problem  studied  here,  the  pressure  under  the  foil  and  the  film 
thickness  are  assumed  to  reach  some  constant  values  p^o  and  h*,  respectively, 
sufficiently  far  away  from  the  bump.  According  to  (1),  then 


_ Tf  - 


f*»  ■ P*  = 

’ ^ 

Neglecting  foil  stiffness,  Eq.  (1)  may  be  rewritten  as 

W4. 


d<*- 


In  the  incompressible  case,  a simple  model  for  the  leakage  factor  is 


U 

(SIDE  LEAKAGE)  ( f “ fj®  j 


(4) 


(5) 


(6) 


where  "h"  is  a characteristic  film  thickness  and  "a"  is  a typical  side  flow  dimen- 
sion (Fig.  2). 


relation 


The  differential  equations  (2),  (5)  must  be  supplemented  by  the 


+ 


h 


(7) 


Assuming  that  the  bump  shape  may  be  approximated  by  a parabola,  w^  is  given 


\ • 4.>|x| 

•»4=( 

Wxi 

J » a-  ' ' ' 

(8) 

0 

where  is  the  bump  protrusion  and  6 its  base  width. 

The  boundary  conditions  for  Eqs.  (2),  (5)  — (8)  are: 
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poo 


It  is  to  be  noted  that  if  the  level  of  h*  is  prescribed  at  -«>,  one  is  not  free  to 
specify  it  at  x -►  »»  (due  to  the  side  flow). 

NORMALIZATION 

The  following  normalizations  will  be  used 


s=t 


If  we  neglect  time  dependent  inertia  effects  as  seen  from  the 
moving  coordinate  ( ^^"C>in  Eqn.  (5)),  and  assume  incompressible  flow,  w( 
find  using  the  above 


where  0 is  a proportionality  leakage  parameter 


and  e,  the  speed  parameter 


& = 


t»/<  ( Ut  - 0,.')  g 
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The  bump  shape  is 


W 


r “ 


as>  lil 

os<  l%l 


The  boundary  conditions  for  steady  state  are; 

P-pP  p'Vo 

4 > <M>  P V o F ' V o 


The  steady  state  solution  of  the  problem  may,  thus,  be  expressed  as 


P = fit  6,W 


(151 


(16) 


RESULTS  AND  DISCUSSION 

Computer  programs  have  been  developed  to  study  the  model  des- 
cribed above.  In  this  report  a limited  range  of  results  of  interest  are  presented. 
Further  results  will  be  published  in  future  reports.  The  present  results  are  based 
on  the  following  assumptions. 

1.  Planar  flow 

2.  The  local  side  leakage  is  a constant  times  (p-Poo).  i-e.,  it  is 
assumed  that  H = 1 in  the  side  leakage  term  of  Eqn.  (11). 

3.  Fluid  inertia  effects  are  neglected. 

4.  Only  the  steady  state  foil  inertia  effects  (in  terms  of  the 
moving  coordinate)  are  included. 

5.  Foil  stiffness  is  neglected. 
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6.  Fluid  compressibility  is  neglected. 

7.  The  head  shape  is  a parabola  of  specified  height,  intersecting 
the  drum  with  a specified  width  6. 

Two  sequences  of  Figures  represent  the  main  results.  Figures  (4a)-(4c)  show 
steady  state  foil  contours  in  relation  to  head  shape  under  various  conditions.  Figures  (6)-(8) 
summarize  the  effects  of  various  parameters  on  the  minimum  head  to  tape  separations. 

In  Figure  {4a),  the  effect  of  e is  studied.  For  concreteness,  one 
may  think  of  e as  a speed  parameter.  As  the  speed  increases,  the  film  thickness 
in  the  rear  of  the  head  grows,  whereas  in  front  of  the  head,  the  film  thickness 
decreases  due  to  the  "inertial  impact"  of  the  tape,  or,  putting  it  differently,  due 
to  the  tape  "getting  in  the  way"  of  the  moving  head. 

It  is  to  be  noted,  however,  that  as  the  value  of  U|-U|.  goes  through 
the  propagation  speed  of  transverse  waves  in  the  foil  j ^ , (referred  to  loosely 
as  "sonic  speed"),  the  sign  of  e changes  (Figure  5).  As  long  as  the  head  shape 
is  symmetrical,  the  solution  for  a given  +e  is  a reflection  about  the  head  axis 
of  the  solution  for  -e  (see  Fig  (4a)  for  e = i 100).  This  fact  clearly  hints  to 
useful  variations  achievable  by  non-symmetric  heads. 

In  Figures  {4b)  and  7,  the  effects  of  side  leakage  are  demonstrated. 

By  increasing  /3,  the  side  leakage  parameter,  the  film  thickness  in  the  rear  of  the 
bump  is  reduced.  Thus,  side  flow  provides  a means  for  controlling  head  to  tape 
separation  in  recording.  Finally,  the  effect  of  head  protrusion  over  the  drum  is 
illustrated  in  Figs  (4c)  and  is  summarized  in  Figure  8.  This  factor  is,  again, 
important  in  achieving  the  desired  flying  characteristics  on  recording.  It  is  perhaps 
worth  mentioning  that  for  sufficiently  large  P and  w^pgj^,  head  to  tape  separation 
may  be  reduced  until  contact  is  established. 

One  may  note  that  in  the  "subsonic"  case,  waves  appear  in  front 
of  the  head,  whereas  in  the  supersonic  case  they  appear  in  the  rear.  Figure  9 
shows  the  wavelength  of  the  tape  contour  as  a function  of  P and  e.  It  is  inter- 
esting to  remark  that  for  each  value  of  0,  there  is  a value  of  e below  which  no 
waves  exist. 
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In  conclusion,  this  work  has  resulted  in  the  development  of 
computer  programs  modelling  head  to  tape  interface  in  helical  recorders.  The 
effects  of  speed,  tension,  side  flow  and  head  protrusion  may  be  studied  for 
different  head  geometries  and  for  optimization  purposes. 
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Figure  7 Minimum  Film  Thickness  vs  Side  Leakage 
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Figure  8 Minimum  Film  Thickness  vs  Bump  Height 
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Figure  9 Wavelength  of  Tape  Contour  vs  Speed  Parameter 


Figure  10  Typical  O.D.E.  Solution.  Upper  and  Lower  Bound  on  the 
Function  F and  its  Derivatives  for  Positive  6 
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APPENDIX  - NUMERICAL  TECHNIQUES 


This  appendix  describes  two  numerical  techniques  utilized  for 
obtaining  the  solutions  presented. 

(a)  Steady  State  Asymptotic  O.D.E. 


If  one  limits  interest  to  the  asymptotic  steady  state  behavior,  the 
time  dependent  term  in  Eqn.  (1)  may  be  dropped.  With  the  side  leakage  term 
approximated  as  a — — , the  equation  may  be  integrated  once  to 

‘'4’'  JC 

' 4% 


To  satisfy  the  boundary  conditions  at  ±00  , a linearized  from  of  the  equation 

may  be  utilized,  i.e.. 


+ 


O 


For  , the  solution  of  this  equation  is  used  to  start  numerical  integration 

in  the  - ^ direction,  after  eliminating  the  exponential  components  of  the  solution 
that  do  not  satisfy  the  boundary  conditions.  The  method  of  solution  is  described 
in  greater  detail  in  [11],  where  it  is  applied  to  a different  problem.  Different  trials 
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generate  upper  and  lower  bounds  for  the  solution.  Once  convergence  is  achieved, 
the  upper  and  lower  bound  solutions  deviate  from  one  another  only  for  large 
values  as  shown  in  Figure  10.  For  6<o  , the  solution  is  a reflection  about  ^=0. 

The  main  advantage  of  this  solution  method  is  the  variable  grid  size 
which  is  related  to  the  local  truncation  error,  leading  to  high  accuracy.  Its  dis- 
advantage is  the  sensitivity  of  the  solution  to  low  film  thicknesses  due  to  the 
appearance  of  in  the  denominator. 

(b)  Time  Dependent  Solution 

The  time  history  of  the  foil  contour  may  be  numerically  solved 
from  Eqn.  (11).  As  an  alternative  to  the  technique  described  above,  the  steady  state 
solution  may  be  found  as  the  asymptotic  limit  of  the  time  history.  This  approach 
is  particularly  attractive  since  we  have  at  our  disposal  a precompiler  [12]  which 
permits  the  generation  of  such  a solution  with  relative  ease.  It  should  be  emphasized, 
however,  that  in  the  derivation  of  Eqn.  (11),  the  time  dependent  terms  of  Eqn.  (5) 
have  been  neglected.  Furthermore,  r (Eqn.  10)  contains  the  factor 
which  “reverses"  time  for  "supersonic"  speeds.  Therefore,  Eqn.  (11)  constitutes  a 
true  physical  representation  of  unsteady  phenomena  only  for  low  "subsonic"  speeds. 
For  higher  speeds  and,  especially  for  "supersonic"  speeds  Eqn.  (11)  should  be 
interpreted  merely  as  a mathematical  device  for  obtaining  steady  state  solutions. 

The  present  project  helped  to  focus  on  some  desired  improvements 
in  the  precompiler. 

1.  A three  dimensional  plot  program  has  been  incorporated  into 
the  precompiler  which  helps  assess  how  far  the  behavior  has 
reached  asymptotically.  The  result  of  the  plot  program  is 
illustrated  in  Figure  11. 

2.  The  method  of  discretization  of  the  term 

^ f ^ ^ 

^ significant  effects  on  the  results. 

If  denotes  the  central  difference  operator  of  order  n,  the 
discretization 
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±(h*^')  ^ S'Ch*s^p) 


(AD 


works  well,  whereas  the  discretization 


fH’  ^ 


(A2) 


» 3rtVHf*F  + H*^*F 


introduces  mass  conservation  errors.  A similar  problem  has  been 
reported  in  [13,  14],  The  capability  of  discretizing  according 
to  (A1)  is  being  introduced  into  the  precompiler. 

3.  Finally,  the  need  for  fine  discretization  in  the  fasi  varying 
regions  of  the  problem  vs  coarse  discretization  in  the  slowly 
varying  regions,  is  quite  acute  in  the  present  problem.  Efforts 
are  underway  to  incorporate  in  the  precompiler  a variable 
grid  size  capability. 
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Figure  11  Development  of  Foil  Contour  From  an  Arbitrarily  Prescribed  Initial 
Condition,  as  Seen  From  a Reference  Moving  with  the  Head 


